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Abstract. Strangeness, charmonium and open charm yields in relativistic 
nucleus-nucleus collisions are considered within statistical model approach as 
potential signals of the quark-gluon plasma. 

I. Introduction 

The present status of the quark-qluon plasma (QGP) in nucleus-nucleus 
(A+A) collisions is somewhat uncertain. Experimental data for A+A colli- 
sions with truly heavy beams have become available: Au+Au at 11 A-GeV 
at the BNL AGS and Pb+Pb at 158 A-GeV at the CERN SPS §. A system- 
atic analysis of these data could yield clues to whether a short-lived phase 
with quark and gluon constituents, the QGP, exists during the hot and dense 
stage of these reactions. This question, whether the QGP is already produced 
with the currently operating heavy ion accelerators, is right now vigorously 
debated. There are also convincing hopes to find new reliable evidences of 
the QGP within a few years. The accelerators of a new generation RHIC 
BNL and LHC CERN will start soon to operate. 

Since a long time || statistical models are used to describe hadron mul- 
tiplicities in high energy collisions. Thermal hadron production models have 
been successfully used to fit the data on particle multiplicities in A+A colli- 
sions at the CERN SPS energies (see, e.g. [0, 01 [51]). Due to the large number 
of particles a grand canonical formulation is used for the modeling of high 
energy heavy ion collisions. Recently, an impressive success of the statistical 
model applied to hadron multiplicities in elementary e + + e~, p + p and p + p 
interactions at high energy was also reported |§. However, in the latter case 
the use of a canonical formulation of the model, which assures the exact con- 
servation of the conserved charges, is necessary (see, e.g., and references 
therein) . 
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The temperature parameter which characterizes the available phase space 
for the hadron production is found in these interactions to be 160-190 MeV 
|J. It does not show any significant dependence on the type of reaction 
and on the collision energy (at the SPS energies and higher). This fact 
suggests the possibility to ascribe the observed statistical properties of hadron 
production systematics in elementary and nuclear collisions at high energies 
to the statistical nature of the hadronization process || ||, || . 

II. Strangeness Production 

The enhanced production of strangeness was considered by many authors 
as a potential signal of QGP formation (see, e.g., Ref. ||). The expectation 
was that strangeness production should rapidly increase when the energy 
transition region is crossed from below. The strangeness to pion ratio is in- 
deed observed to increase in A+A collisions. It seems however that this is not 
a signal of the QGP. The low level of strangeness production in p+p interac- 
tions as compared to the strangeness yield in central A+A collisions, called 
strangeness enhancement, can be also understood to a large extent as due 
to the effect of the exact strangeness conservation. The canonical ensemble 
treatment of the strangeness conservation leads to the additional suppression 
factors imposed on the strange hadrons production in small systems created 
in p+p collisions || 0. Another important point is that for the chemical 
freeze-out parameters, temperature T and baryonic chemical potential 
found for the SPS energies the strangeness to entropy ratio is larger in the 
equilibrium HG than in the equilibrium QGP 0. 

To estimate the strangeness to entropy ratio let us consider the quantity 

where N s and N s are the numbers of strange quarks and antiquarks, and 
S is the total entropy of the system. In the QGP we use the ideal gas 
approximation of massless «-, d-(anti) quarks and gluons, strange (anti)quarks 
with m s = 150 MeV and (anti) charm quarks with m c = 1500 MeV. For the 
HG state the values of N s and N§ are calculated as a sum of all s and s inside 
hadrons, and S is the total HG entropy. The behaviour of R s (P for the HG 
and QGP is shown in Fig. 1 as a function of T for hb = 0. In the wide 
range of T = 200 + 500 MeV one finds an almost constant value of R s in the 
QGP which is smaller than the corresponding quantity in the HG. The total 



2 



entropy as well as the total number of strange quarks and antiquarks are 
expected to be conserved approximately during the hadronization of QGP. 
This suggests that the value of R s at the HG chemical freeze-out should be 
close to that in the equilibrium QGP and smaller than in the HG at chemical 
equilibrium. Therefore, the strangeness suppression in the HG would become 
a signal for the formation of QGP at the early stage of A+A collision at the 
CERN SPS energies. 



The same conclusion was obtained in Ref. ||10|| . In the model presented 
in that paper it was assumed that due to the statistical nature of the cre- 
ation process the strangeness in the early stage is already in equilibrium 
and therefore possible secondary processes do not modify its value. As the 
strangeness to entropy ratio is lower in the QGP than in the confined matter, 
the suppression of strangeness production is expected to occur when crossing 
the transition energy range from below. The total strangeness production is 
usually studied using the experimental ratio 

(A) + (K + K) 

Es = — 5o — ' {) 

which measures the ratio of the mean multiplicities of A hyperons and K, K 
mesons to the multiplicity of pions. This is an experimental analog of the 
quantity R s ([!]). In Fig. 2 the experimental data of E s ratio (0) are shown 
together with its theoretically expected behaviour according to Ref. [0 . The 
deconfinement phase transition causes a nonmonotonic behaviour of E s with 
the collision energy. The maximum of E s is in the energy region between the 
AGS and SPS. 

III. C harmonium and Open Charm Production 

Charmonium production in hadronic and nuclear |LlJ collisions is 



usually considered to be composed of three stages: the creation of a cc pair, 
the formation of a bound cc state and the subsequent interaction of this 
cc bound state with the surrounding matter. The first process is calculated 
within perturbative QCD, whereas modeling of non-perturbative dynamics is 
needed to describe the last two stages (see, e.g., |H| and references therein). 
The interaction of the bound cc state with matter causes suppression of the 
finally observed charmonium yield relative to the initially created number of 
bound cc states. This initial number is assumed to be proportional to the 
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number of Drell-Yan lepton pairs, which then allows for the experimental 



study of the charmonium suppression pattern. It was proposed |1J] that 
the magnitude of the measured suppression in nuclear collisions can be used 
as a probe of the state of high density matter created at the early stage of 
the collision. The rapid increase of the suppression (anomalous suppression) 



observed when going from peripheral to central Pb+Pb collisions [15] is often 



attributed to the formation of the QGP ||16 |. 



It was recently found [|T(], [L7| that the mean multiplicity of J/ip mesons 



increases proportionally to the mean multiplicity of pions when p+p, p+A 
and A+A collisions at CERN SPS energies are considered. We illustrate this 
unexpected experimental fact by reproducing in Fig. 3 the ratio (J/ip)/ (h~) 
as a function of the mean number of nucleons participating in the interaction 
for inelastic nuclear collisions at the CERN SPS. The (J/ip) and (h~) denote 
here the mean multiplicities of J/ip mesons and negatively charged hadrons 
(more than 90% are ir~ mesons), respectively. 

In the standard picture of the J/ip production based on the hard cre- 
ation of cc pairs and the subsequent suppression of the bound cc states the 
observed scaling behavior of the J/ip multiplicity appears to be due to an 
'accidental' cancelation of several large effects. A very different picture was 



proposed in Ref. [18| which explains a scaling property of the J/ij) multiplicity 



(h-) 



const(A) (3) 



by assuming that a dominant fraction of J/ip mesons is produced directly at 
hadronization according to the available hadronic phase space. J/ijj mesons 
are neutral and unflavored, i.e., all charges conserved in the strong inter- 
action (electric charge, baryon number, strangeness and charm) are equal 
to zero for this particle. Therefore, its production is not influenced by the 
conservation laws of quantum numbers. Consequently, the J/ip production 
can be calculated in the grand canonical approximation and, therefore, its 
multiplicity is proportional to the volume, V, of the matter at hadronization. 
Thus, the statistical yield of J/ip mesons at hadronization is given by 

W) = (23+ J V i^P exp|(p * 1/2/rg| _ j (4) 



^^g), (2j + 1 ),(^f 



exp , . 



4 



where j '• = 1 and = 3097 MeV are the spin and the mass of the J /if) 
meson and Th is the hadronization temperature. The previously mentioned 
results of the analysis of hadron yield systematics in elementary and nuclear 
collisions within the statistical approach indicate that the hadronization tem- 
perature Th is the same for different colliding systems and collision energies. 
This reflects the universal feature of the hadronization process. 

The total entropy of the produced matter is proportional to its volume. 
As most of the entropy in the final state is carried by pions, the pion multi- 
plicity is also expected to be proportional to the volume of the hadronizing 
matter. Thus the scaling property ([!]) follows directly from the hypothesis of 
statistical production of J /if) mesons at hadronization and the universality of 
the parameter Th- Since elements of hadronizing matter move in the overall 
center of mass system the volume V in Eq. (^J) characterizes in fact the sum 
of the proper volumes of all elements in the collision event. 

The hypothesis of statistical production of J /if) mesons at a constant 
hadronization temperature Th leads to the prediction of a second scaling 
property of the J /if) multiplicity, namely: 

S const(^S) (5) 

which should be valid for sufficiently large cm. energies, y/s. This scaling 
property is illustrated in Fig. 4 which shows the ratio (J /if)) / (h~) as a func- 
tion of y/s for proton-nucleon interactions. The experimental data on J/ip 
yields are taken from a compilation given in The values of (h~) are 



calculated using a parameterization of the experimental results as proposed 

Onwards from the CERN SPS energies, y/s = 20 GeV, the ratio (J/if))/{h~ 
is approximately constant, in line with the expected scaling behavior The 
rapid increase of the ratio with collision energy observed below y/s = 20 GeV 
should be attributed to a significantly larger energy threshold for the J /if) 
production than for the pion production. In terms of the statistical approach 
the effect of strict energy-momentum conservation has to be taken into ac- 
count by use of the microcanonical formulation of the model. 

The statistical J /if) multiplicity @ depends on two parameters, Th and 
V. However, a simple way to estimate of the crucial temperature parameter 
in Eq. (||) from the experimental data is possible, provided that we find a 
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second hadron which has the properties of the J /if} meson, i.e., it is neutral, 
unflavored and stable with respect to strong decays. The best candidate is 
the rj meson. The multiplicity of 77 mesons seems to obey also the scaling 
properties @ and fl5|). The independence of the (77) / (n°) ratio on the collision 
energy was observed quite a long time ago pl| . Recent data on 77 production 



suggest that (rj) j (tt°) ratio is also independent of the size of the colliding 
objects. In central Pb+Pb collisions at 158 A GeV the ratio (t7)/(vt ) 



0.081 ± 0.014 is measured [p2j| . It is consistent with the values of the ratio 
reported for all inelastic p+p at 400 GeV H (0.077±0.005) and S+S at 
200 A -GeV J24| (0.12±0.04). From the measured ratios, {J/^)/{h~) and 
(77) / (7r°) , we estimate a mean ratio (J/if>)/(r)) = (1.3 ± 0.3) • 10~ 5 . Here 
we use the experimental ratio (tt°) / (h~) = 1 in N+N interactions P?5| ]. 
Under the hypothesis of the statistical production of J /if) and 77 mesons at 
hadronization the measured ratio can be compared to the ratio calculated 
using Eq. (f|): 

W) ^ 3 m\ K 2 (m^/T H ) 
(77) m2 K 2 (mjT H ) ' 1 ' 

where m v = 547 MeV is the mass of the 77 meson. This leads to an estimate of 
the hadronization temperature, Tjj ~ 176 MeV. A graphical solution of Eq. 
(|6[) is shown in Fig. 5 which illustrates the high sensitivity of the estimate 
of the Th parameter by using the (J/ip)/(r)) ratio. This is due to the large 
difference between mass of the J /if) and the 77 mesons as the right hand side 
of Eq. (|) is approximately proportional to exp[(m ?? — m^)/T H ]. 

In the proposed model the creation of the J/if> mesons is due to the 
straight thermal production at hadronization and not due to the coalescence 
of cc quarks produced before hadronization. Therefore the yield of J /if) 
mesons is independent of the production of open charm, which is carried 
mainly by the D mesons in the final state. The D mesons multiplicity is 
determined by the number of cc quark pairs created in the early parton stage 
before the hadronization. Assuming chemical equilibrium of charm in the 
QGP stage of A+A collision we can estimate the ratio of the open charm 
hadrons to pions. The charm to entropy ratio is defined by the quantity 

similar to Eq. ([!]) for the ratio of strangeness to entropy. The behaviour of 
R c (0) for the HG and QGP is shown in Fig. 6 as a function of T for fig — 0. 



6 



In Fig. 6 the quantity R s ([I]) from Fig. 1 is also presented for a comparison. 
The behaviour of R c is completely different from that of R s . R c in the QGP 
is strongly increasing with T and its values are much larger than in the HG. 
The experimental analog of R c is E c = (D)/(tt), which measures the ratio 
of the mean multiplicities of D mesons to the multiplicity of pions. The 
assumption of the conserved total entropy and the total number of charm 
quarks and antiquarks during the hadronization of QGP leads then to the 
picture of chemical non-equilibrium hadron gas at the chemical freeze-out 
with a strong enhancement of D mesons yield. The ratio of D mesons to 
pions should strongly increase with the collision energy from the SPS to the 
RHIC. 

IV. Summary 

Statistical production of starngeness and charm discussed in this talk can 
be summarized as the following: 

• The transition to the QGP in A+A collisions could be seen as a non- 
monotonic dependence on collision energy of the strangeness to pion 
ratio in the energy region between the AGS and the SPS. 

• The yield of J ftp mesons at the CERN SPS can be understood assum- 
ing the statistical production of J/ip at the hadronization, and it is 
sensitive to the hadronization temperature. The D mesons multiplic- 
ity is determined by the number of cc quark pairs created in the early 
QGP stage. The ratio of D mesons to pions is expected to increase 
strongly with the collision energy from the SPS to the RHIC. 
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Figure 1: R s (1) at /ib = for the HG (solid line) and the QGP (dashed 
line) . 
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Figure 2: The dependence of the strangeness to pion ratio, E s (2), for 
the central A+A collisions (closed circles) and nucleon-nucleon interactions 
(open squares) as a function of the collision energy measured by the variable 

F = (y/s — 2m7v) 3//4 /v / ^ 1 ^ 4 - The solid line shows predictions of the statistical 
model of Ref. [10|. A transition to the QGP takes place between the AGS 
(F w 2) and the SPS (F w 4) energies. 
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Figure 3: The ratio of the mean multiplicities of J/i[) mesons and negatively 
charged hadrons for inelastic nucleon-nucleon (square) and inelastic O+Cu, 
O+U, S+U and Pb+Pb (circles) interactions at 158 A-GeV plotted as a 
function of the mean number of participant nucleons. The results for Pb+Pb 
interactions measured for different centralities of collisions are shown by open 
circles. For clarity the N+N point is shifted from (N P ) = 2 to (N P ) = 5. 
The dashed line indicates the mean value of the ratio. 
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Figure 4: The ratio of the mean multiplicities of J ftp mesons and negatively- 
charged hadrons for inelastic proton-nucleon interactions as a function of the 
collision energy in the center of mass system. 
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Figure 5: The (J/ip)/ (rj) ratio calculated under hypothesis of the statistical 
production of J/ip and r\ mesons at hadronization (solid line) as a function 
of the hadronization temperature. Band shown by dashed lines is drawn at 
±<t around the mean experimental value of the (J/ip)/(rj) ratio. The dotted 
line indicates T H = 17Q MeV. 
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Figure 6: R c (0) at ji B = for the HG (solid line) and the QGP (dashed line). 
R s behaviour from Fig. 1 is also presented (upper lines) for a comparison. 
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